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We have computed correlation functions for nucleons and extracted the masses for positive- and negative- 
parities. Use of group theory plays an important role in obtaining sources that have good overlap to higher 
spin states and minimum contamination from unwanted states. In the simulation three distinct sources and 
corresponding sinks that transform according to the Gi irreducible representations are tested and used to form 
matrices of correlation functions. Diagonalizations give us clear mass splittings between low-lying states and 
excited states for both parities. 



1. INTRODUCTION 

Lattice QCD calculations of baryon masses are 
based on the analysis of 2-point correlation func- 
tions, defined by 



Ci i (t,to) = 5^(0|B i (x,t)fl J -(0,t )|0> 



(1) 



where Bi are color-singlet baryon source fields. 
Summation over x projects to zero total momen- 
tum. For large t the contribution from the lowest 
mass eigenstate dominates. An effective mass is 
calculated from M e g = In I „?^ t '* ^ . 



hedral group" , or 2 0. The reduction of the group 
2 falls into three spinorial IRs: G±,G2, and H 
with respective dimensions 2, 2, and 4. Tabled 
lists the inverse relation of the reduction of SU(2) 
to 2 pp. Fermionic operators on a cubic lattice 



Table 1 

Correlation of IRs and total angular momentum 



\C{t+x,t ) / 

Ideal sources maximally overlap the states of 
interest: (n| Bi(to) |0) ~ S n i, where n is an inter- 
mediate state. By using the orthogonality rela- 
tion of basis vectors spanning irreducible repre- 
sentations (IRs) of the cubic group, one can at- 
tempt to realize this relation by taking Bi to be a 
linear combination of source operators that have 
a good overlap with state n. 

There are 48 discrete spatial rotations that con- 
serve the cubic symmetry of spin 1 /2 objects, and 
they form a group, called "double covered octa- 



states. 


IRs 


3 


Gi 


1/2,7/2,9/2,11/2,... 


G 2 


5/2,7/2,11/2,... 


H 


3/2,5/2,7/2,9/2,... 
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can be written in terms of basis operators that 
span IRs of 2 0. Since lattice calculations pre- 
serve cubic symmetry, different IRs do not admix. 
However, each IR is an admixture of different to- 
tal angular momenta. Total angular momentum 
is not a "good" quantum number on the lattice 
because continuous rotational symmetry is bro- 
ken. It is expected that higher spin states have 



2 



higher masses, therefore within a given IR the 
lowest spin state dominates at large t. 

2. NONLOCAL SOURCES 

A general expression of three-quark baryon op- 
erator would be written as 

B M (x) = h,h,h, M, a, & 7 )e o6c 

Ur'q a a(h)U b /q b ;(l 2 )U^q^(h), (2) 

where ot,/3, and 7 are Dirac indices, a, b, c, a', 
b', and c' are color indices, and li , I2 , and I3 are 
spatial displacement vectors. ' a is the gauge 
link connecting the spatial path from x to li. The 
coefficients tie together a set of spatial displace- 
ment vectors L, I2, and I3, and Dirac indices a, (3, 
and 7 in such a way as to yield gauge invariant 
operators that transform according to IRs of 2 
and definite isospin. 

In order to have sources that are commensu- 
rate with the size of baryons, it is appropriate to 
smear the sources, using for example, Gaussian 
smearing, 
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(3) 



where a is the smearing width, N is an integer, 
and D 2 is the gauge-covariant three-dimensional 
Laplacian operator. For different radial distribu- 
tions, additional powers of Laplacian operators 
are used, e.g., F 2 = D 2 F , F4 = D 2 F 2 , and so 
on. 

The simplest nonlocal operator has L = L> = x, 
I3 = x + aei. Applying spatial rotations of the 
cubic group to this operator leads to 6 x 8 = 48 di- 
mensional representations, where six is for the di- 
rections of aei and eight for the three spin degrees 
of freedom. Reduction of these representations 
results in 6 Gi's, 2 G2's, and 8 H's. From these 
"one-link" constructions, we used both positive- 
and negative-parity nucleon sources from G\ ba- 
sis in the simulation. 

More complicated sources have been con- 
structed by choosing the initial operator as, for 
instance, li = L> = x, I3 = x + aei + ae 2 in 
Eq. ©. 



3. SIMULATION DETAILS AND PRE- 
LIMINARY RESULTS 

For a test of these sources, we formed two 3x3 
matrices of correlation functions with positive- 
parity and negative-parity, by using projection 
operators 1 ± 74. All operators transform ac- 
cording to Gi IRs. Isospin is selected to be 
I = I z = 1/2. The operators used in the sim- 
ulation are shown below, 

s[f 2 = [d(C l5 )u] (1 ± 74) (1 + 7375 K (4) 



D z s 



2j±) 



1/2 = [d(C l5 )u] (1 ± 74 )(1 + ~f 3l5 )D 2 u, (5) 

p[f 2 = [d{C lb )u] (1 T 74)(1 + 7375)£ + w 

+ KC75M (1 T 74)(1 - l3l 5 )D z u, (6) 

where Di is a covariant derivative and D + = 
D x + iD y . Color and Dirac indices are implicit 
here. The upper sign is for positive-parity and 
the lower sign is for negative-parity. All quark 
fields are smeared using Eq. ©. Equation (@J 
is the widely used G\ irreducible source, which 
mainly couples to the ground state of the nu- 
cleon: we call it "s^" even though from Ta- 
ble ^one sees that there can be admixtures with 
spin 7/2 or higher. Equation JSJ is s-wave with 
an additional Laplacian operator acting on the 
third quark. Equation JHJ is constructed from 
one-link sources as described above. The deriva- 
tive operators and spin-projection operators in 
Eq. 10 that act on the third quark can be rewrit- 
ten as -D+u_!/ 2 + D z Ui/2- This corresponds to 
the Clebsch-Gordan formula for a p 1 / 2 state, i.e., 



P1/2 = -\l g^llX-l/2 



2^10X1/21 



(7) 



where Y\ m represents a discretized version of 
a spherical harmonic, and x±\/2 represent up 
and down two-component spinors. Because of 
this correspondence, Eq. © is called the "p^" 
source. 

We use 16 3 x 32 isotropic lattice with periodic 
boundary conditions. The smearing width that 
is used in Eq. © is a = 4.5, with N = 50. 
We use 221 configurations 2 generated from the 



2 Configurations are from the NERSC lattice archive. 
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isotropic Wilson gauge action at (3 = 6.0 and the 
Wilson fermion action with k — 0.1550, which 
corresponds to pion mass about 542 MeV. 
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Figure 2. Effective masses of negative-parity Gi 
states. 



Figure 1. Effective masses of positive-parity G\ 
states, related to the diagonal correlators men- 
tioned in the text. 



We adopt a variational method to extract ef- 
fective masses |3I4I5| . Correlation matrices are 
diagonalized at a time slice t v . This produces 
three orthogonal eigenvectors, each of the form 
v = (vi,V2,v 3 ), where Uj's are the coefficients of 
Syl, D 2 s^j\, and Py\, respectively. The result- 
ing eigenvectors are then used to calculate "di- 
agonal" correlators of the form v T C(i)v at other 
times. For positive-parity the eigenvectors ob- 
tained at t v — 10 give the most clear mass split- 
ting, which is shown in Fig. (JJJ. Filled diamonds 
have v « (0.44,0.57,0.68), i.e., each source is 
almost equally weighted. This state reaches a 
plateau for 5 < t < 10, suggesting that the 
ground state of the nucleon has been obtained. 
Filled circles and open boxes are for other orthog- 
onal eigenvectors. The effective masses based on 
the other orthogonal eigenvectors are larger than 
the ground state mass, indicating that excited 
states of positive-parity are being resolved. 



For negative-parity, the eigenvectors at t v = 
6 give the most clear mass splitting, which is 
shown in Fig. J5J). Filled diamonds have v = 
(0.37, 0.56, 0.71), where the contribution from the 
p[/\ state is significant. The effective mass of this 
state is about 40% larger than the positive-parity 
ground state; the corresponding iV*(l/2 - ) phys- 
ical state is about 63% heavier than the nucleon. 
Filled circles and open boxes have larger error 
bars, but they indicate larger masses than dia- 
monds. 

The present results are preliminary. Work is 
under way to improve the chiral properties of 
the action and to use a large set of sources with 
anisotropic lattices. 
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